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1 Introduction 



Most of interesting physical models and theories are described by the gauge invariant Lagrangians, which are 
singular and in Hamiltonian formulation lead to the constrained systems [1]. For the constrained Hamiltonian 
systems there are in principle two waysj^ of quantization [1-2]: 

1 . "First quantize and then reduce" . 

2. "First reduce and then quantize". 

The present paper deals mainly with the reduction procedure of the latter. For gauge theories this procedure 
is a restriction on the constraint surface Al and then farther reduction to the physical phase space M = M./G, 
which is the space of gauge orbits. 

If the action of the gauge group (G) on the constraint surface [M) is regular, then the manifold of orbits 
[M = M-jG ) is well defined and it possesses the symplectic structure. Coordinates on M. are gauge invariant, 
true physical degrees of freedom. 

Quite often, in practical applications this theoretical scheme of reduction encounters technical problems 
related to the explicit construction of M. with its symplectic structure: apart from the mathematical difficulties, 
the physical content of the true degrees of freedom may be quite unpredictable. 

Widely used practical reduction scheme is a gauge fixing procedure by some conditions g) = [1],[3]. 
For simple cases the explicit form of the true physical variables is obvious and this reduction scheme works 
perfectly. But in general, as it was shown in [4] (namely for the Yang-Mills theory), the space of gauge orbits 
{M. /G) cannot be obtained by "simple" gauge fixing. Problem of gauge fixing, of course, reflects the above 
mentioned possible non-trivial structure of a physical phase space [1] , [5] . 

Another reduction scheme can be based on the gauge invariant variables (GIVs) [1-2]. As a rule, the GIV 
is constructed from the structure of gauge transformations. If one can find the complete set of GIVs, then it 
allows to describe the physical phase space M-jG with its symplectic structure. This paper deals with such 
gauge invariant approach by using of restricted 1-forms. We also consider situations when only a part of GIVs 
is known. Analysis of a structure of restricted 1-form helps to find the remaining part of GIVs. 

Note that the reduction scheme with 1-forms for arbitrary constrained systems was proposed in [6] . In these 
papers elimination of extra variables was based on the Darboux theorem. Sometimes Darboux theorem is not 
effective in applications and choice of GIVs is just a practical way for the realization of this reduction program 
for gauge theories. 

The paper is organized as follows: In Section 2 the reduction scheme to the GIVs is introduced and, for the 
illustration, simple examples are considered. One more example of (2 + l)-dimensional massive photodynamics 
is given in Appendix. In Section 3 this scheme is applied to the finite dimensional system with SU (2) gauge 
group of symmetry. This system can be considered as a toy model of the Yang-Mills theory with fermions. It 
is shown that there is an essential difference between this SU{2) and the corresponding U{1) model. Structure 
of GIVs in SU{2) case can be interpreted as the confinement phenomenon. In Section 4 we study infinite 
dimensional model, where gauge group is any semi-simple one. The GIVs are constructed and full reduction is 
accomplished. It is shown that the model is equivalent to the Yang-Mills theory with some boundary conditions. 
Final section is for remarks and conclusions. 



2 Reduction scheme in gauge invariant variables 

Starting from the gauge invariant Lagrangian L = L[qk,qk) {k — l,...,iV) and using the Dirac's procedure 
[la], or the first order formalism [6] we arrive to the extended phase space F with coordinates {pk,<lk) and the 

^In this paper wo do not consider the path integral approach 
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action 

S = j Pkdqk - [H{p, q) + XaMP, q)]dt, (2.1) 
k = l,...,N; a = l,...,M; {N > M), 

where (j)a{p,Q) are constraints, H{p,q) is a canonical Hamiltonian and are Lagrange multipliers. The con- 
straint surface — is defined by 

Mp,Q)=0 (2.2) 

and the following relations are fulfilled: 

{H, </.Jr = d%, {cPa, <Pb}r = flb^c- (2.3) 

Index r on the left hand side indicates that Poisson brackets are calculated on the extended phase space. 
Function ^ = £,{p, q) is called GIV [1] if and 

{^,<f>a}r = dt<f>b, (2.4) 

where \m denotes restriction on M and functions rf^ (as well as and /^^ in (2.3)) are assumed to be regular 
in the neighbourhood of A4 . 

Each GIV — ^ has the class {^} of equivalent GIVs on F [1]. A gauge invariant function ^ is equivalent 
to ^ if = . On the other hand, the function is a constant along the gauge orbit (on {A4)) and it 
defines the function f on the physical space M = A4 /G. Thus (_\m and f denote the GIV — ^ in different 
context. If there is no ambiguity, we will use the notation ^ for all of them. 

Maximal number of GIVs (2.4), which are functionally independent on the constraint surface A4, is 2{N—M) 
[lb]. Suppose that : a = 1, 2(A'' — M)} is such complete set of GIVs. Then one can prove [6] that 

1. Pkdqk\M = 6*1 + O2, with 

a) dOi = 0, 

h) 02 = 0a{OdC, (2.5) 
c) detujap 7^ 0, where oJafiiC) = da^iS — dpQa'i 

2. H{p,q)\M=h{0- 

Main statement of (2.5) is that after restriction on the constraint surface M., dependence on extra (non-physical) 
variables is present only in the term 9i , which is a "total derivative" . 

Since d9i = 0, it gives no contribution to the variation of a restricted action. We can neglect it and for the 
reduced system get 

S\m = S = J 0a{OdC - Ki)dt. (2.6) 
Hence dynamics for GIVs is described by the Hamilton equations 

e=a;«^(05/3/i(0. (2.7) 

where w"''(^) is the inverse to the symplectic matrix ci;„/3 = daO^ — d^ffa and it defines the Poisson brackets of 
the reduced system 

{C,f}M=oj'-'^HO- (2.8) 

So the reduced system (2.6)-(2.8) is an ordinary (non-constrained) Hamiltonian system which can be quantized. 
It should be noticed that, in general, any 2(A'' — M) number of GIVs are only local coordinates on the 

physical phase space Ai and respectively (2.5)-(2.8) have the local meaning. Global description can be achieved 
by the set of GIVs which defines the global structure of physical phase space M. Number of such GIVs is 
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greater than 2{N — M), but on the constraint surface there are relations among them. Just these relations 
define the geometry of M. For the illustration let us consider the following example of (2.1)-(2.3) [7a]: 

S = j P- dq- + MMdt. (2.9) 

Here p and q are vectors of R^, canonical Hamiltonian is zero, 

(j)i=p-q, (j)2=p'^q'^ - [p- qf -r'^ 

and r is a parameter. These constraints are Abelian ({0i, = 0) and the second constraint (/>2 can be written 
in the form 

where J = (f X p is an angular momentum. 

It is clear that the physical phase space is two dimensional and the components of angular momentum J 
are GIVs (they commute with constraints, since constraints are 0(3) scalars). On the constraint surface these 
three components arc related by J - J = and define the physical phase space Ai as the two dimensional sphere. 
So any two coordinates (as well as the 1-forms 9x and O^) are only local ones (on the phase space geometry of 
constrained systems see [5a]). 

Described reduction scheme ((2.5)-(2.8)) can be used when all 2{N — M) GIVs are known. For the practical 
application of the scheme one can introduce any additional (to GIVs) variables r]^,..., r]^ to have the coordinate 
system 

(^\...,e^(^-^),r?\...,r?^) 

on A4. Then calculating restricted 1-form Pkdqk\M in these coordinates and taking its differential wc can find 
the symplectic form w = i^a/3(^)d^" A d£^^ . Note that in practical calculations it is possible to single out the 
1-form 6*2 = OaiOd^" and arrive to (2.6). 

Application of this procedure to the model (2.9) gives 02 = zd^p, where z and <p are the cylindrical coordinates 
on the sphere: 

Ji = \J — z'^cosif, J2 = \J — z'^sirup, J3 = z. 

It is clear that zd^p is not the global 1-form, but its differential has a continuation to the well-defined symplectic 
form on the sphere [8] 

^ _ Ji{dJ2 /\dJ3) + Mdh A dJi) + J:i{dJi A dJ^) 

After this the system can be quantized by the geometric quantization [9] (see also [7a], [10b, c]). Consistent 
quantum theory exists only for the discrete values of the parameter r. 

Generalization of the scheme to the infinite dimensional case is straightforward (in Appendix we present the 
example of massive photodynamics in (2 + 1) dimensions). If we use the Dirac's observables [11]: 

= e^^"'(^^)t/; (2.10) 

in the ordinary QED, we will easily obtain photons in the Coulomb gauge and the "four-fermion interaction" 
for the "dressed fermions" (compare to the example in Section 3 and see [6], [12]). 

Note that the commutation relations of the complete set of GIVs (2.8) can be derived by calculations of 
Poisson brackets on the extended phase space too [1]. This more standard procedure is based on the fact that 
the Poisson brackets of any two GIVs is again GIV. Indicated procedure and the scheme described in this 
paper ((2.5)-(2.8)) arc almost equivalent. Only, sometimes, calculations of differential forms is more practical 
(especially, when the canonical quantization is not applicable [9]). 

In general, from the structure of gauge transformations one can easily find only some part of GIVs and 
construction of the complete set (2.5) is troublesome. In many cases, our approach with differential forms, can 
be effectively used for the solution of this problem too. 

Let us consider situation when we know the set of GIVs {^'^ : /i = 1, iiT}, where N — M < k < 2{N — M). 
We can introduce any additional variables r]^, 7y2JV-M-if have the coordinate system on A4 and calculate 
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the restricted 1-form Pkdqk\M- Suppose that we can single out "total derivatives" and the differentials d£,^ in 
the form 

Pkdqk\M = dF{^,rj) + e^{^,r])d^^'. (2.11) 

Then using (2.5) wc can easily conclude that Of^{(_,ri) will be GIVs. Note that passing to the GIVs — 
helps to get the form (2.11). For the illustration of this method we apply it to the relativistic particle [Ic], with 
the 1-form 6 = pdq — podqo and the constraint surface M: p^ — w? = 0, (po > 0). The momenta p are gauge 
invariant and after restriction on A4 we have 

= pdq - Vp"^ + m^dqo. 
One can easily rewrite it in the following form 

0\m = d{p- \/p^l^qn) - {q j=^==qo)dp. 

Evidently the coefficients of the differentials —dp axe GIVs, canonically conjugate to p: 

a/P ^ + 

Gauge invariance of Q also can be established from the relation 

L = Vp 2 + m?Q, (2.12) 

where L arc the generators of Lorentz transformations. Since all generators of the Poincarc group (P^jM^^) 
are GIVs, the same property have the coordinates Q. On the constraint surface p^ — m? = 0, (po > 0) all these 
are the functions only of the reduced variables p and Q. 

Reduced system can be easily quantized in momentum representation: p = p and Q = ihV. Operator 
ordering problem arises only for the generators (2.12). Then the standard Lorentz covariant measure of a scalar 
product 

<*2|*1 >= /^J^=*2(p)*l(p) 

J ^p^ + ' 



corresponds to the ordering L = ih^p"^ + m^V. 



3 The finite dimensional models with and SU{2) gauge symme- 
tries 



In this section we consider the finite dimensional model with SU(2) gauge group of symmetry. From the 
beginning it is difficult to see all GIVs and we use the method described at the end of Section 2. Obtained 
structure of GIVs is quite unexpected. For comparison we present the corresponding U{1) model too. These 
U(l) and SU{2) models can be considered as the toy models of the electrodynamics and the Yang-Mills theory 
(with matter), respectively. In classical description all "fields" are assumed to be c- numbers. 



A. The model with U{1) symmetry 



Let us consider the action 



S = I dt[Uiptp - i^i)) - mi^tp + Aoitp^J - kE) +EA- ]-E\ (3.1) 
J 2, z 
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where all "fields" {■tp, tp, Aq, A, E) are functions only of the time t; m and k {k ^ 0) are parameters. Similarity to 
the electrodynamics is apparent from the notations. At the same time (3.1) has the form (2.1), where Aq = \{t) 
is a Lagrange multiplier and ^ = "iptp — kE is a constraint (we use the time derivatives instead of differential 
form when it is convenient). 

Non-zero Poisson brackets are 

{i^A}=i, {E,A} = 1 

and we get the gauge transformations 

V'(t) — > e+^"WV(i). ^(*) — ' e-'"^*^{t), 

A{t) — y A{t) + ka{t), E{t) — y E{t). (3.2) 



Then 



Ao{t) — (*) + "(*) 



leaves the action (3.1) invariant. 

Reduced system is two dimensional and two GIVs can be chosen as 



*inv = e^""^' *inv = e*^^^ (3-3) 



(compare to (2.10)). Here, the reduction procedure (2.5) is trivial and we get 



S= dt 



2 (^inv^inv ^inv^inv) ™*inv*inv (*inv^inv) 



(3.4) 



So the "gauge field" A has vanished and physical excitations are only the "dressed fields" ^'i^v (with "four- 

fermion interaction"). 

This model has a simple generalization in case of multi-component gauge field A with gauge transformations 

A — > A + ka, 

where k are parameters (fc^ ^ 0). The GIV ^'jny constructed similarly to (3.3) (or (2.10)). Then, after 

reduction, "longitudinal" (to the k) component of the gauge field A vanishes and physical variables are only 
"transverse" ones and the constructed "dressed field" "^[ji^y So, for these Abelian models, the structure of 
GIVs is very similar to the physical observables of the electrodynamics [6], [12]. 



B. The model with SU{2) symmetry 



For the model with SU (2) gauge group of symmetry we consider the action 



-I 



dt 



^{'ipa-lpa - V'aV'a) - m^'aV'a + Ao{j + J) + EA - ^E"^ 



(3.5) 



Here tpa are 2-componcnt spinors (a = 1, 2), m is a parameter, A and E are three dimensional vectors, Aq 
lagrange multipliers and the angular momenta j and J are given by 



j = i)-tp, J = AxE, 



(3.6) 



where a are the standard Pauli matrixes. 

Connection with the Yang-Mills theory is obvious. 
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Non-zero Poisson brackets are 

{■4'a,'4'l3} = iSa0, {Em,An} = Smn, (m,n) = 1,2,3 (3.7) 

and the constraints ^ = j — J generate the gauge transformations: 

ujip, tp ipio~^, A — > ujAu)~^, E — > wEw"-^, 

where oj{t) G SU{2) and 

1 -. 1 - 

A = -Aa, E = -Ea. (3.8) 

Then for Aq = ^A^a wc get Ao uAqw^"'^ — ilulo^^ . 

Any scalar product of the vectors A,E, J,j will be GIV. But on the constraint surface {j + J = 0) only 
three of them are functionally independent. 

If we choose these independent GIVs as: 

l„ = l{A^ + E% h = l{EA), l^ = l{A^-E% (3.9) 
then from (3.7) we get the SL{2,R) algebra: 

{IM = ^i^^-^flp, where gi'" = diag {+,-,-,). (3.10) 

Since there are three constraints, the physical phase space is 4-dimensional. To construct the fourth GIV 
and find the full symplcctic structure we use the method of Section 2 (see (2.11)). 

For the parameterization of the constraint surface we introduce new variables (j, /i, </)): 

i = ^(/M+/i2), /i = ^(/ii - /l2), 



where 

Then for the 1-form we get 



$ = V-l + (/92, <^ = </'l-V2, (3.11) 

V'a = \/^e-''^", i>a = VKe'''^ (a = 1,2). 



% — — 

-{ipadtpa - V'a#a) = jd^ + hdifi. (3.12) 



The vector j (3.6) in these new coordinates will take the form 



y/j^ — h'^COSif 

y/ p — h? sirup 
h 

and = 2^. Note that on the constraint surface we have (see (3.9)): 1*^1^, = and for the fixed j the 
commutation relations (3.10) define well known symplectic structure on this hyperboloid (see e.g. [13]). 
If we introduce the ortho-normal basis (cj • = ^ik, x e*, = e^fcCfc ) : 




ei = I COS(p I , 62 = -y 



costp 
sirup 



J 

, , , 63 = T, 

h 



then A and E can be parameterized as follows: 

A = eiQi + 62^2, E = eipi + 62^2, 
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where 

Piq2-P2qi=j- (3.13) 
Calculating the restricted 1-form EdA\_M in these new coordinates and using (3.13) we obtain 

EdA\M = pidqi +P2dq2 - hd(p. (3-14) 

Comparing (3.12) and (3.14), we sec that there is a cancellation of the 1-form hdip. This means that the 
corresponding degrees of freedom vanish. 

Now, it is convenient to introduce the polar coordinates for two-vectors (91,32) and {pi,P2)- 

qi = rcosP, pi = pcos'y, 

q2 = rsinP, p2 = psiwj. 

Then three of them (r, p and (/? — 7)) are connected with the GIVs (3.9): 

r"^ = 2{la + l2) = l+, p"^ =2{lo-l2) = l-, rpcos{l3 - j) = 2h. 
Using these relations we finally get the reduced 1-form 

e\j^=jd^ + h^ where 1? = $-/?. (3.15) 

'+ 

So the coordinate 1? = $ — ,9 is the fourth GIV. Respectively the reduced Hamiltonian takes the form 

H\m = 2mj + 7 ^ (3.16) 

'+ 

and this is the complete reduction. 

Note that the second part of the reduced 1-form Zid(lnZ_|_) defines the above mentioned symplectic structure 
on the hypcrboloid /^Z^ = \,f [13]. 

We see that the physical picture of this reduced system essentially differs from the corresponding Abclian 
case. Here, after reduction some part of degrees of freedom of the "gauge field" {A), as well as some part of 
"matter field" (tp) degrees of freedom have vanished. Below we shall see that in quantum theory vanishing of 
"matter field" degrees of freedom can be interpreted as the confinement phenomenon. 

Geometric quantization [9] is a natural way for the construction of quantum theory of the reduced system 
(3.15-3.16), but in principle one can use canonical quantization too. For this aim it is convenient to introduce 
(global) "creation" and "annihilation" variables 

a+ = ^je'^, a = yTje-^^ (3.I7) 

and in quantum theory we get the discrete values of j = a^a. Then, quantization of the system with the 
canonical 1-form Zid(ln/+) and the Hamiltonian (3.16) (for obtained discrete values of j), gives the irreducible 
representations of 51/(2, R) group (see e.g. [13a]). 

Next, from (3.11) we have the relation N = 'ipai'a = 2j. It is natural to interpret the corresponding operator 
(N = 2j) as the ^ particle number operator. In quantum theory we have 

[N, a+] = 2a+ 

where a+ is a physical creation operator (3.17). So in physical excitations (created by the operator a+) there 
arc states only with even number of "fermions" . This fact also can be seen from the structure of the variable 
a+ (see (3.17) and (3.11)). It has the phase factor e^^'^^~^'^^\ So in quantum case it will create (see [14]) the 
pairs of "dressed" V'-particles. 

Note that for the similar finite dimensional constrained systems such "confinement" -like phenomenon has 
been derived by the "first quantize and then reduce" method (see [15]). In that approach reduction of the 
extended "Hilbert" space by the conditions ^aj^p/iys >= forbids the states with certain quantum numbers. 
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4 Field theory models with non-Abelian gauge group of symmetries 



For the finite dimensional models of the previous section the gauge group G acts on the configuration space of 
"gauge field" {A) and on the phase space of "matter field" {ij)). This is the standard situation for Yang- Mills 
theories. 

Using notations (3.8) we have 

EdA=<'E.,dA>, (4.1) 

where < , > is a scalar product in corresponding Lie algebra A. Thus, the Lie algebra A can be interpreted 
as the configuration space of a "gauge field" A and trivial cotangent bundle as the phase space. 

If one takes a manifold of semi-simple Lie group (G) as the configuration space, then there are the natural 
actions (left and right) of G on this manifold and one can construct similarly the gauge theory where phase 
space is the cotangent bundle [16] T*G = {{g, R)\g & G,R G A}. On T*G the symplectic form is given by 

uj = de, with e =< R, g~^dg > . (4.2) 

Generators of the left and right transformations (g — > ojg, g — > gui) are respectively left and right currents 
(L = gRg^^, R). Choosing gauge transformations as the right action, we get that constraints are (p = R = 0. 
So the "gauge field" part in the action takes the form 



/ 



< R,g~'^dg >-{<A,R> +H{R,g))dt, (4.3) 



where A S ^ is a Lagrange multiplier, and H some gauge invariant Hamiltonian. 

Field theory generalization of (3.5) is the standard Yang-Mills theory. In this section we consider corre- 
sponding generalization of (4.3) with the action 



S = J dt J d°-^x^Yl <Rk,gk^gk>+e<Ao,4>>^ -H 



(4.4) 



where gk{x,t) G G and Rk{x,t),Ao{x,t) £ A; H in gauge invariant Hamiltonian, Aq — Lagrange multipliers, 
4>{x,t) = eX^^/ Rk{x, t) constraints, e — coupling constant (see below). 

The "1-form" J2k=i Rk,9k^dgk > defines the equal time Poisson brackets (see e.g. [16]): 

{Rk,ai^)'^i,b{y)} = Ski6{x-y}f^^Rk,c{x) 
{9k{x),Ri,aiy)} = 5ki6{x - y}{gkTa{x)) 

{9k{x),gi{y)} = (4.5) 

where the set {Ta\ Ta G A} forms any basis of Lie algebra. Ra =< Ta, R > and the last two relations are matrix 
equalities [16]. So for the constraints (f>a =< Ta, (f) > wc have 

{MS!), My)} = 6{x-y)ra,M^). (4.6) 
Corresponding gauge transformations are 

9k — > 5feW, Rk — > w~'^RkW (4.7) 
and one can easily construct GIVs such as 

9ki = gkgi^ and Lfe = gkRk9k^- (4-8) 
The Hamiltonian H in (4.4) is any functional of such GIVs. 
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Since (4.8) gives sufficient number of GIVs we can use the scheme described in Section 2. The first non- 
trivial case is the 3-dimensional space-time. If we introduce g = gig^^ as the variables and Ri,R2 and g2 
as the 77 variables of the scheme (see (2.11)), then for the "1-form" 6 =< Ri, g^^dgi > -|- < R2, g2^dg2 > 
(integration over is assumed) we immediately get 

< Ri +R2,g2^dg2 > + < g2Rig2^,g~^dg > 

and after reduction we have 

e\M=<r,g-^dg>, (4.9) 

where r = g2Rig2^ is also GIV. 

So the structure of the 1-form is the same, only the number of variables was reduced. One can check that 
this is true for other dimensions too. 

It is clear that the phase spaces of the systems with 1-forms (4.1) and (4.2) are essentially different and they 
cannot be transformed to each other. But in field theory when one has the infinite number of such spaces there 
is a non-local transformation (see [17c]): 

Ak^^g^^dugk, Ek = -egi:'^d^:^{Lk)gk, (4.10) 

which transforms the system (4.4) into the Yang-Mills theory with the same gauge group G. Indeed, from (4.4) 
and (4.10) one can check that 

D-1 

^i) = ^ dkEk + e[Ak, Ek] (Gauss law) 
fc=i 

and 

< Ek,Ak >=< Rk,gk^gk > + (total derivatives) (4-11) 
To get the corresponding Hamiltonian of the Yang- Mills theory [17]: 

^ . D-1 ^ D-1 

H=^ d'^-^xiY, < Ek,Ek > +2 H < Fki,Fki >) 

fe=l k,l=l 

with Fki = dkAi — diAk + e[Ak, Ai], one has to choose in (4.4) 

= d^'-'ie^ < d^'Lk,d^'Lk >+^< dk{gkidigik),dk{gkidigik) >]■ (4.12) 

So one can assume that the system (4.4) with the Hamiltonian (4.12) is equivalent to the ordinary Yang-Mills 
theory with some boundary conditions (which allow to invert (4.10) and neglect the total derivatives in (4.11). 

Boundary behaviour is a subtle problem even for simple models of field theory (see e.g. the Appendix). For 
the Yang-Mills theory it is too complicated and we do not consider it here. 

Unfortunately the complicated form of the Hamiltonian (4.12) does not simplify after the reduction proce- 
dure. For example, for the considered 3-dimensional case the reduced Hamiltonian takes the form 



H =^ j d^x < V > < 82^1,82^1 > + 



4 < di{gd2g-^),diigd2g-' >], (4.13) 

where I = grg^^ . 

Gribov's ambiguity problem has stimulated many papers on the gauge invariant description of the Yang- 
Mills theory and the reduced system (4. 9), (4. 13) is the one possible version (literature and new results on this 
problem see in [le]). The main problem of such approaches is a complicated form of the Poincare generators in 
GIVs [17]. For example, the Hamiltonian (4.13) is non-local in fields an non-analytical in coupling constant. 
So the standard perturbative quantization is not applicable. 
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Note that such Hamiltonian with corresponding symplectic form was obtained in [17c] by Dirac's brackets 
formalism. 



5 Conclusions 



Of course, there was an essential progress in the study of constraint systems since the paper [18], but from the 
point of view of practical applications still there is no universal approach. The method presented in this paper 
is a one possible practical procedure towards the quantization of gauge theories. 

As it was mentioned in the introduction, for the gauge invariant systems there is an alternative way of 
quantization when one "first quantizes and then reduces" . In general, there are two problems in such approach: 

a) construction of physical states l^phys > as the solutions of the equations (f)a\^phvs >= 0, where ^ are 
constraint operators. 

b) problem of scalar product for the physical states. 

Sometimes the first problem is only a technical one (for the Yang-Mills theory see [18]), but in general both 
this problems are related and need further investigation [20]. 

In this paper we have not mentioned other important methods such as the path integral approach [lb], [3] 
and BRST quantization [21] (for a review see [10a]). 

Quantization procedure is not unique even for the ordinary, non-constrained systems [8] , [22] . It depends on 
the choice of canonical variables (if they globally exist), operator ordering, etc. Therefore it is not surprising that 
different quantization procedures of constraint systems generally lead to the non-equivalent quantum systems 
[3],[23]. 

As it was mentioned in Section 2, for a reduced classical system generally there are no global canonical co- 
ordinates and usual canonical quantization is not applicable. This, together with technical problems of classical 
reduction, was the main obstacle in general formulation of the approach "first reduce and then quantize" . 

Geometric quantization [9] and other "new" quantization schemes [7], [10a], [24] allow to quantize Hamiltonian 
systems without global canonical structure too. At the same time essential progress was made in the construction 
of classical reduction schemes [6] . Therefore for the wide class of constrained systems the quantization method 
"first reduce and then quantize" seems to us to be technically preferable. Here it should be mentioned about 
the possible combination of the two quantization schemes: if a reduced classical system is complicated, then on 
the cotangent bundle of a reduced phase space one can construct new extended system with simple constraints 
and next use the first way of quantization [7],[10a,b]. Of course, the question, which is the "correct" quantum 
description of a given classical system, remains open. 
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Appendix 

The 2 + 1 dimensional massive photodynamics is described by the Lagrangian (see e.g. [25]): 

1 777 

C^--F^,F'^- ^-e'^'^^F^^A^. (A.l) 
We choose 17^,^ = diag{+, — , — ), e°^^ = 1 and in the first order formalism [6] obtain 

J dtj^^ d^x[{E, - ^e,jAj)A, - ^{E,E, + B^) + Ao{d,E, - mB)], (A.2) 

where 

E, = Fh = A,~ 9,^0, B = -e„F,, (e,, = e"^') 

and we have neglected the boundary term iPxdi\Aii{^€ijAj — Ei)\. 

If we use "1-forms" instead of time derivatives (see comment after (3.1)), then the action (A.2) takes the 
form (2.1) with Aq playing the role of Lagrange multiplier. 

For the reduction we choose E\ and E2 as the variables A\ as the additional variable r\ (see (2.11)) and 

get 

-E^E^ - \\E,E, + \(duEkf\ + 
^m, 2 777^ dt 

where 

e = i[SiAi + S2X(Ai + l£;2)] 

2 777 

and the operator K = d^^d2 is assumed to be symmetric due to the corresponding boundary conditions. 

Neglecting the term as the total derivative, we get the local Hamiltonian theory with the canonical 
commutation relations 

{i;2(x), Si (y)} = 777(5(2) (x- 7/) (A 4) 

and the quadratic Hamiltonian 

i/ d^x[E,E, + \{dkEkf]. (A.5) 

2 777^ 



dt I d^x 

fl2 



(A.3) 
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The energy- momentum tensor also can be expressed only through Ei and E2: 

Too = llEiEi + ^{dkEkf], Toi = -CijEjidkEk). (A.6) 

Let us briefly stop on the boundary conditions. We can assume that a boundary behaviour of the physical 
variables {Ei,E2) should provide the Poincare invariance of the reduced system (A.3)-(A.6), while a boundary 
behaviour of the fields of the initial system (A.l) should allow the outlined reduction procedure. 

Generators of the Poincare group (constructed from the energy-momentum tensor (A.6)) generate transfor- 
mations of El and E2 according to the Poisson brackets (A. 4). The class of functions Ei{x) and E2{x) should 
be invariant under these transformations. It is natural to choose the class of smooth, rapidly vanishing at the 
infinity functions. 

For the diagonalization of the Hamiltonian and momentum let us take the Fourier transformation: 

Ej{x)=iJd'p^^Ej{p) (A.7) 
and introduce the longitudinal and transverse components: 



Ej{p) = ^^ei{p)-'-^e2{p), (A.5 



where \p\ = v^pfTpf- 

Then diagonalization of the energy and momentum will be achieved in the variables 



a{p) = 



^ei(rt +ze2(p) ^_,^(^^ 
^ei(-rt-ze2(-p) ^.^(p) 



a*{p) = " ^ ' ^^_^e"^W^ (A.9) 



'2ojp 



with 



bl 



Note that for the chosen class of Ei{x) and E2{x) the longitudinal and transverse components of Ej{p) have 
the singularity at the origin [p = 0) and the phase factor e^'^'^^-' is necessary to cancel it. On the other hand 
one can easily check that the class of smooth functions a(p), a*{p) is Poincare invariant. Just this phase factor 
was introduced in [25] to avoid anomalies in the commutation relations of the Poincare algebra of quantum 
operators. As we have seen this phase factor is connected to the Poincare invariance of the classical system too. 

After description of the class of physical variables one can go back and find the class of gauge potentials Afj,. 
One can show that these classes for massive and ordinary photodynamics in (2+1) dimensions are different. 
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